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Abstract 

We provide a significant extension of the Hyperboloidal Foliation Method introduced by the authors in 2014 in 
order to establish global existence results for systems of quasilinear wave equations posed on a curved space, when 
wave equations and Klein-Gordon equations are coupled. This method is based on a (3+1) foliation (of the interior 
of a future light cone in Minkowski spacetime) by spacelike and asymptotically hyperboloidal hypersurfaces. In the 
new formulation of the method, we succeed to cover wave-Klein-Gordon systems containing “strong interaction” 
terms at the level of the metric, and then generalize our method in order to esstablish a new existence theory for the 
Einstein equations of general relativity. Following pioneering work by Lindblad and Rodnianski on the Einstein 
equations in wave coordinates, we establish the nonlinear stability of Minkowski spacetime for self-gravitating 
massive scalar fields. 


Stabilite nonlineaire globale de l’espace-temps de Minkowski pour les champs massifs. 

Resume 

Nous generalisons la Methode du Feuiletage Hyperboloidal introduite par les auteurs en 2014 pour traiter des 
systemes quasilineaires couplant des equations d’ondes et des equations de Klein-Gordon. Dans cette nouvelle 
formulation, nous reussissons a traiter des termes metriques “d’interaction forte”. Nous appliquons cette methode 
pour demontrer la stabilite nonlineaire de l’espace de Minkowski pour les equations d’Einstein des champs scalaires 
massifs auto-gravitants. En suivant un travail de Lindblad et de Rodnianski, nous analysons la structure des 
equations d’Einstein en coordonnees d’ondes, qui constituent precisement un systeme d’equations d’ondes quasi¬ 
lineaires avec “interaction forte”. 
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Version franqaise abregee 


Nous etudions le probleme de l’existence globale en temps de solutions regulieres d’equations d’ondes 
nonlineaires, avec deux objectifs principaux : 

- Nous generalisons la methode du feuiletage hyperboloidal [8] proposee par les auteurs en 2014. 

- Cette methode nous permet d’analyser les equations d’Einstein et de demontrer la stability nonlineaire 
de l’espace de Minkowski en presence d’un champ scalaire massif auto-gravitant. 

Les aspects nouveaux de notre methode (voir le texte en anglais pour plus de details) sont les suivants : 

(1) inegalites de Sobolev et de Hardy associees au feuilletage hyperboloidal de l’interieur d’un cone de 
lumiere (de l’espace-temps de Minkowski). 

(2) Estimations de type pour l’equation d’onde et pour l’equation de Klein-Gordon sur un 

espace courbe. 

(3) Une hierarchie d’estimations d’energie ayant des croissances algebriques en temps differentes. 

Par ailleurs, en ce qui concerne le traitement des equations d’Einstein proprement dit, nous combinons 
differentes idees et techniques, enoncees comnre suit : 

(4) decomposition de la courbure de Ricci (“null forms”, “quasi-null forms”), 

(5) decomposition de tenseurs basee sur la condition d’onde, 

(6) integration le long de caracteristiques, et 

(7) une hierarchie d’energies adaptees a la structure des equations d’Einstein-Klein-Gordon. 


English version 
1. Introduction 

We are interested in the global-in-time existence problem for small amplitude solutions to nonlinear 
wave equations posed on a curved spacetime, with a two-fold objective. First, we provide a significant 
extension of the Hyperboloidal Foliation Method, proposed by the authors [8] in 2014. This method is 
based on a 3 + 1 foliation (of the interior of a future light cone in Minkowski spacetime) by spacelike 
hyperboloidal hypersurfaces and on Sobolev and Hardy-type inequalities adapted to this hyperboloidal 
foliation. This method applies to a broad class of coupled nonlinear wave-Klein-Gordon systems on curved 
space, and takes its root in a work by Klainerman [7] and Hormander [4] on the Klein-Gordon equation. 

In comparison to our earlier formulation [8], we are now able to encompass a broader class of coupled 
systems involving “strong interaction” terms (as we call them, see below). Recall that Klainerman used 
the decomposition of the wave operator along an hyperboloid al foliation and was able to establish a 
(non-sharp) decay rate of f -5 / 4 and conclude with a theory of existence of quasilinear Klein-Gordon 
equations. Later on, Hormander worked directly with the energy on hyperboloidal hypersurfaces and 
Sobolev inequalities and derived the sharp rate f -3 / 2 for the problem posed in flat space. 

Our work [8] provides a way to combine both approaches and encompass systems of coupled wave and 
Klein-Gordon equations. Importantly, we work directly within the hyperboloidal foliation and, in order 
to encompass equations posed on curved space, we establish sup-norm bounds leading us to the sharp 
rate f -3 / 2 . 

Our second objective is to apply this method to the Einstein equations of general relativity and, by 
constructing spacelike and asymptotically hyperboloidal hypersurfaces, we offer a new strategy of proof 
in order to establish the nonlinear stability of Minkowski spacetime. Our method applies self-gravitating 
massive scalar fields, while all earlier works were restricted to vacuum spacetimes or spacetimes with 
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massless scalar fields; cf. the pioneering work by Clrristodoulou and Klainernran [3], the proof by Lindblad 
and Rodnianski [13,14] (based on the wave gauge), and the extension by Bieri and Zipser [1], 

One of the simplest wave-Klein-Gordon model is, in flat space, Qit = P(du,dv), + v = Q(8u,dv), 
where P, Q are quadratic forms in the first-order derivatives du = (d a u) and dv = (d a v) and the two 
unknowns u,v are defined over Minkowski space R 3+1 . (Here a = 0,1,2, 3.) Many models arising in 
mathematical physics involve interactions between massive and massless fields. Let us mention the Dirac- 
Klein-Gordon equations, the Proca equation (massive spin-1 field in Minkowski spacetime), the Einstein- 
massive field system, and the field equations of modified gravity described by the Hilbert-Einstein func¬ 
tional f(R g ) dv g with, typically, f(R g ) = R g + k(P s ) 2 and k > 0, where R g is the scalar curvature of 
a Lorentzian manifold ( M,g ). 

The vector field method was introduced by Klainernran [6,7] around 1980-1985. It primarily applies to 
quasilinear wave equations posed on the (3 + l)-dimensional Minkowski spacetime, and leads to global-in¬ 
time well-posedness results when the initial data are sufficiently small in some Sobolev spaces. The method 
relies on the use of the conformal Killing fields of Minkowski space, suitably weighted energy estimates, and 
the so-called Klainerman-Sobolev inequalities. Nonlinearities are assumed to satisfy the ‘null condition’ 
and a bootstrap argument is formulated and relies on time decay estimates. In comparison, quasilinear 
Klein-Gordon equations have atracted less attention in the literature, despite pioneering contributions by 
Klainernran [7], Shatah [15], and Hormander [4]. 

In our work [8,9,10,11], we have addressed this major challenge of developing a method for quasilinear 
wave-Klein-Gordon systems on curved space. The main difficulty conres from the fact that a smaller 
symmetry group is available to deal with Klein-Gordon equations, since the scaling field tdt + rd r does 
not commute with the linear Klein-Gordon operator in flat spacetime. While additional decay for Klein- 
Gordon equations, that is, t~ 3 ' 2 in four dimensions is available (solutions to wave equations decaying only 
like f _1 ), a robust technique to deal with the coupling of wave equations and Klein-Gordon equations 
is required and we have developped the Hyperboloidal Foliation Method precisely for that purpose. For 
earlier contributions on the analysis of Klein-Gordon equation with a limited number of Killing fields, we 
refer to Katayama [5] and the references cited therein. 


2. The hyperboloidal foliation method 

Following [8], we rely solely on the Lorentz boosts (or hyperbolic rotations), which generate a foliation 
(of the interior of a light cone) of Minkowski spacetime by hyperboloidal hypersurfaces (that is, surfaces 
of constant distance from a base point). We introduce here Lorentz-invariant energy norms based on 
these boosts and are able to revisite all the standard arguments of the so-called vector field method and 
carefully analyze the energy flux on the hyperboloids. With a suitable extensions of Sobolev and Hardy 
inequalities adapted to the hyperboloidal foliation, we are able to encompass a broad class of coupled 
systems. 

The method in [8] is based on the hyperboloidal hypersurfaces 1~L S := {{t,x) j t > 0; t 2 — \x\ 2 = s 2 } 
parametrized by their hyperbolic radius s > so > 1: and consider the foliation of the future light cone 
K, : = { (t, x) / |a:| < t— l}. Note in passing that s < t =% s 2 . We imposed initial data on the hypersurface 
t = So > 1 or directly on the hyperboloid s = sq. Our energy estimates therein are formulated in domains 
limited by two hyperboloids, that is, /C[ S0)Sl ] := {(t, x) /\x\ < t - 1, (s 0 ) 2 t 2 - \x\ 2 (si) 2 ,f >0}. Our 

analysis is performed in the semi-hyperboloidal frame (as we propose to call it), consisting of the Lorentz 
boosts L a := x a dt + td a , a = 1,2,3 and a time-like vector. More precisely, by definition, this frame 
consists of the following three vectors tangent to the hyperboloids d a := and the timelike vector 
d 0 := dt- Accordingly, we have the semi-hyperboloidal decomposition of the (flat) wave operator: Qu = 
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— yrd 0 8 0 u — — ^-d a d 0 u + 2 a —a— — | dtU. (In comparison, the standard choice in the literature 

is the ‘null frame’, containing three vectors tangent to the light cone.) 

The hyperboloidal energy associated with the hypersurface TL S involves certain weighted derivatives, 
and we want to point out that we will use the full expression of the corresponding energy flux on the 
hyperboloids. Let us also mention one important functional inequality. 

Lemma 2.1 (Sobolev estimate on hyperboloids) For all functions u defined on a hyperboloid 1~L S in 
Minkowski space R 3+1 and with sufficiently fast decay, one has sup( tx ) eWa t 3 ? 2 \u(t, x)| < 2 |/|ss2 ||£ 7 u!| l 2 (h.) 
(for s ^ so > 1) with summation over L e {L a = x a dt + td a }, where I denotes a multi-index. 

The hyperboloidal foliation method is based on a hierarchy of bounds for the curved metric and the 
source-terms. Specifically, in [8], we used three levels of regularity and algebraic growth rates and, re¬ 
markably, our bound is uniform for the low-order energy of wave components. 


3. The Einstein-massive field system 

We present a new method for proving the nonlinear stability of Minkowski spacetime, which applies 
to self-gravitating massive scalar fields [10,11]. The statement of the problem is as follows (following 
Choquet-Bruhat et al. [2]): we search for a spacetime ( M,g ) satisfying the Einstein equations 

-ftotfl 2 ^ ct/3 

for the stress-energy tensor of a scalar field f>, that is, 

T a p : = VafiVpfi - + V(<j)))g a p, 

where the potential is taken to be V((j>) := ^-</> 2 (c > 0 being the mass of the scalar field). Using the 
contracted Bianchi identities, it is not difficult to derive the Klein-Gordon equation \Z\ g <t> = = c 2 0. 

Our objective is to study the associated Cauchy problem when the initial data set is a perturbation of a 
spacelike lrypersurface in Minkowski space. 

Theorem 3.1 (Nonlinear stability of Minkowski space for massive fields) Consider the Einstein- 
scalar field system in wave coordinates (that is, \Z\ g x a = 0): 

\Jg9a0 = (Qap + P a p) {g\ dg,dg) - 2(d a (j)dp<t> + V(4>)g a p), n g </> - V'(</>) = 0, 

where D g := g a & d a ’dp’ is the so-called reduced wave operator, Q a p are null terms, and P a p are “weak 
null” terms. Consider an initial data set ( M,~g,k ) which is close to a spacelike slice of Minkowski space 
and is asymptotically hyperboloidal (in a suitable sense) Then, the initial value problem for the Einstein- 
massive field system admits a global solution in wave coordinates, which defines a future geodesically 
complete spacetime (M, g ). 

In [10,11], the massive field is assumed to be compactly supported, but this assumption is not essential 
and is removed in [12], which also treats the more general theory of modified gravity. Our proof relies on 
the wave gauge, after the pioneering work of Lindblad and Rodnianski [13,14]. The following challenges 
and techniques are in order: 
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Tensorial structure: The geometric structure of the Einstein equations in combination with the wave 
coordinate condition □ g x a = 0 allows us to decompose the quadratic nonlinearities as a sum of null 
terms and “weak null” terms. 

- Hyperboloidal foliation: Having fewer Killing fields at our disposal, we rely on the foliation generated 
by the Lorentz boosts, that is, the hyperboloids of Minkowski space and we introduce Lorentz- 
invariant energy norms. As explained in Section 2, we need to establish Sobolev and Hardy-type 
inequalities on hyperboloids. 

Sharp pointwise estimates: We derive estimates for, both, the wave equation and Klein- 

Gordon equations on curved space. We use a technique of integration along well-chosen curves (see 
below). 

- Hierarchy of energy bounds: Several levels of regularity and time growth rates are required in our 
bootstrap argument, and successive improvements of the estimates are performed in the proof. 

In the rest of this text, we present some of our technique for a simplified, but challenging, model. 


4. Wave-Klein-Gordon model with strong interactions 

Consider the following system with strong interactions at the metric level, which (we formally extract 
from the Einstein equations in wave coordinates): 

—□« = P a Pd a vdpv + Rv 2 , —\Jv + uH al3 d a dpv + c 2 v=0. (1) 

Theorem 4.1 (Nonlinear wave-Klein-Gordon model with strong interaction) Consider the non¬ 
linear wave-Klein-Gordon model (1) with given constants P Q/3 , R, H a $ and c > 0. For any N > 8, there 
exists e = e(N) > 0 such that if the initial data satisfy ||(itOi K))||.fl'- W + 1 (R3) + ||(wi, Ui)||#iV( R 3) < e then the 
Cauchy problem for (1) admits a global-in-time solution. 

We proceed with a bootstrap argument based on a hierarchy of energy bounds posed along the hyper- 
boloidal foliation (stated in (2), below). The proof of Proposition 4.2 relies in particular on a suitable 
decomposition of the Klein Gordon equation on curved space as well as an involved ODE lemma. 

Proposition 4.2 (T 00 —L°° estimate for Klein-Gordon equations on curved space) Consider the 
Klein-Gordon equation on a curved background —\Z\ g v + c 2 v = f with metric g a @ = m Q/3 — h a ^ given by 
a perturbation of the Minkowski metric, and with compactly supported data prescribed on a hyperboloid 
v\u SQ = vq, dtv\u SQ = vi for sufficiently smooth and spatially compactly supported data Vo,v\. Then, in 
the future of TL So , one has 

s 3 ^ 2 \v(t,x)\ + - s 3 ^ 2 \dj_v(t, a:)| < V(t,x) 
s 

with V defined below and := dt + ^j-d a . 

Note that d ± is orthogonal to the hyperboloids for the Minkowski metric and coincides, up to an 
essential factor 1/f, with the scaling vector field. We use the notation ht, x ( A) := h (At/s, Xx/s) (with 
s 2 = t 2 — r 2 ) and consider the derivative in A, that is 

h' t ( A) = -dth°°(Xt/s,Xx/s) + - — d a h°°(Xt/s,Xx/s) = - d ± h°°(Xt/s,Xx/s). 

’ s s s 

Fix a constant C > 0 (chosen later on) and define the function V first “far” from the light cone 0 < r/t < 
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V(t,x):= f||uo||i“(«» 0 ) + ll u ilU°°(H S0 )) (l + f \K,x(s)\e C ^ lht ’ xW]dX da) 

Jso 

+F(s) + f F(s)\h'JA)\e c SIKJ»l^ds 

Js 0 

s 2 — l 

and then “near” the light cone < r/t < 1 by 

V(t,x) := F(s) + f F{s)\ti t x (s)\e c ^K*W\ dx ds 
Js(r/t) 

with S(r/t) := J Here, F is dehned by a suitable integration of the given source-term f. 


Proposition 4.3 (T 00 —L 00 estimate for the wave equation with source) Let u be a spatially com¬ 
pactly supported solution to the wave equation — Qu = / with vanishing initial data and source f satisfying 
I/I 5; (t 3 s 2/ for some exponents 0 < p < 1/2 and 0 < \v\ < 1/2: 


r 1 1 


\u{t,x) I < { 


up, (t — rY M t’ 

1 (t-rY 

w t i+v ’ 


0 < v 1/2, 
-1/2 ^ v < 0. 


The proof is based on the explicit representation formula for the wave equation. We now sketch our 
bootstrap argument, based on the following assumptions (with | J\ = k ): : 


E(s, d I L J u) 1 / 2 + s~F 2 E{s, d I L J v ) 1 / 2 Cyes™, \I\ + \J\ s: N, 

E{s, d I L J u) 1/2 < Cie, \I\ + \ J\ < IV - 4, (2) 

E(s, d z L J u) 1/2 < Cies kS , \I\ + | J\ < N - 4. 

Using Sobolev and Hardy inequalities adapted to the hyperboloids, we then deduce basic decay bounds. 
Next, we derive the following sharp sup-norm bounds for |/| + | J| < TV — 4, which are at the heart of our 
argument: 

supt|L J u| + sup {t/s) ll2 ~ AS t :il2 \d I L J v\ + sup (t/sY /2 ~ 45 t 3/2 \d ± d I L J v\ < Cy£s u . (3) 

Hs Us Us 

We proceed as follows: 

- First bound for the wave component (L 00 -!! 00 estimate for wave equations) for |/| + | J\ < N — 7: 
\d I L J u\ < C-yct-* 12 + (C 1 £) 2 {t/sY (k+4)S t~ 1 s (k+4)s . 

- Second bound for the wave component and first bound for the Klein-Gordon component 
for wave and K-G equations)) 

\u(t,x)\ < Cyet~ l \ |u| + -\dj_v(t,x)\ < Ci£(t/s)~ 2+7S s~ 3 ^ 2 . 

s 

- Second bound for the Klein-Gordon component (again the L 00 -!! 00 for K-G) for |/| < N — 4: 

|d x d 7 u(t,a')| < C'ie(t/s) _ 3 / 2 + 4 l 5 t _3 / 2 ; |d 7 u(t,2;)| < Gie(t/s)- 1 / 2 + 45 r 3/2 . 

Third (and sharp, except for the higher order derivatives of v) bound for the wave and Klein-Gordon 
components for |/| + | J| ^ TV — 4: 

sup (f|L J 'u|) < Ci£S kS , 
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sup((Vs) 3_75 s 3/2 |d_ L <9 / L' 7 u|) + sup Ut/s) 2 ~ 75 s 3 ^ 2 \d I L J v\) < C l£ s kS , 

Hs Hs 

sup Ut/s) 1 - 7S s 3 / 2 \d a d I L J v\) < C!es k5 . 

H a 

The pointwise estimates follow from Propositions 4.2 and 4.3, the bootstrap assumptions, the structure 
of the equations, and various commutation properties enjoyed by the vector fields under consideration 
in the Hyperboloidal Foliaiton Method. Finally,we can conclude and close our bootstrap argument by 
returning to the (differentiated) system 

-Ud I L J u = d T L J (P afi d a vdpv) + 8 1 L J (. Rv 2 ) , 

+ uH afi d I L J v + c 2 d I L J v = ~[d I L J ,uH aP d a dp]v, 
and showing that all source-terms provide integrable contributions to the energy. 
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